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ABSTRACT 

Data collected by the Pierre Auger Observatory (Auger) provide evidence for anisotropy in the arrival directions 
of cosmic rays (CRs) with energies >57 EeV that suggests a correlation with the positions of active galactic 
nuclei (AGN) located within ~75 Mpc and 3.2° of the arrival directions. This analysis, however, does not take 
into account AGN morphology. A detailed study of the sample of AGN whose positions are located within 
3.2° of the CR events (and extending our analysis out to ^150 Mpc) shows that most of them are classified 
as Seyfert 2 and low-ionization nuclear emission-line region (LINER) galaxies whose properties do not differ 
substantially from other local AGN of the same types. Therefore, if the production of the highest energy CRs 
is persistent in nature, i.e., operates in a single object on long (> Myr) timescales, the claimed correlation 
between the CR events observed by Auger and local active galaxies should be considered as resulting from a 
chance coincidence. Additionally, most of the selected sources do not show significant jet activity, and hence, 
in most conservative scenarios, there are no reasons for expecting them to accelerate CRs up to the highest 
energies, ~10 20 eV. If the extragalactic magnetic fields and the sources of these CRs are coupled with matter, 
it is possible that the deflection angle is larger than expected in the case of a uniform source distribution due to 
effectively larger fields. A future analysis has to take into account AGN morphology and may yield a correlation 
with a larger deflection angle and/or more distant sources. We further argue that the nearby radio galaxy NGC 
5128 (Cen A) alone could be associated with at least 4 events due to its large radio extent, and PKS 1343-60 
(Cen B), another nearby radio galaxy, can be associated with more than 1 event due to its proximity to the 
Galactic plane and, correspondingly, the stronger Galactic magnetic field the ultra high energy CRs (UHECRs) 
encounter during propagation to the Earth. If the UHECRs associated with these events are indeed accelerated 
by Cen A and Cen B, their deflection angles may provide information on the structure of the magnetic field in 
the direction of these putative sources. Future 7-ray observations (by, e.g., Fermi Gamma-ray Space Telescope 
[Fermi] formerly Gamma-Ray Large Area Space Telescope [GLAST], and High Energy Stereoscopic System 
[HESS] in the Southern hemisphere) may provide additional clues to the nature of the accelerators of the 
UHECRs in the local Universe. 

Subject headings: cosmic rays — galaxies: active — galaxies: individual (Cen A, Cen B, PKS 2158-380, 
PKS 2201+044) — intergalactic medium 



1. INTRODUCTION 

The spectrum, origin, and composition of CRs at the highest 
energies (> 10 18 eV = 1 EeV; hereafter UHECRs) has been 
a puzzle since their discovery almost 50 year s ago (for a re- 
view see lCroninl2005tlN"agano & Watsonl2000h . The isotropy 
of the arrival directions of UHECRs above 10 18 eV suggests 
an extragalactic origin, though Galactic objects such as fast 
rotating neutron stars with ultrastrong magnetic fields ("mag- 
netars" ~) are capable of accelerating particles up to ^10 20 eV 
dHillaslll984l) . The UHECR energy losses due to photopion 
production on the cosmi c microwave background (CMB) - 
the so -called GZK effect (iGreisenll 1 9661 IZatsepin & Kuz ' mini 
1966) - mean that the sources of the highest energy particles 
should be cosmologically close, within ^100 Mpc. However, 
observations of the spectrum at these high energies are ex- 

1 Also Hansen Experimental Physics Laboratory, Stanford University, 
Stanford, CA 94305 

2 Also Astronomical Observatory, Jagiellonian University, ul. Orla 171, 
30-244 Krakow, Poland 



tremely challenging due to the low overall event rate (~1 per 
km 2 per century at ~10 20 eV). Fine structur e in the CR spec- 
trum above 10 18 eV has been predicted (e.g jBerezinsky et al.l 
2006) - the GZK cutoff, a pile-up bump, and a dip due to pho- 
toproduction of pairs - and indeed found in experiments, but 
the interpretation is not straightforward due to the large un- 
certainty in the source distribution, their injection spectrum, 
and CR chemical composition. 

Quite a few extragalactic objects discussed in the lit- 
erature are potentially capable of accelerating UHECRs. 
Among those are shocks from the epoch of the large scale 
structure formation, 7-ray bursts, galaxy clusters, AGN 
(in particular, th ose AGN with powerful jets) and power- 
ful AGN flares (iFarrar & Gruzinovl 120081) . the lobes of gi- 
ant radio galaxies (Stanev 2007a) , and new ly born magne- 
tars in other galaxies (Ghisellini et al.1 120081) : various ex- 
otic top-down s cenarios have also been disc ussed but seem 
unlikely (e.g., lAbraham etail l2007al l2008bl) . Extragalac- 
tic jets and their extended radio lobes have been proposed 
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as o ne of the most likely acceleration sites of UH ECRs 
(e.g. iBiermann & Sttittmatterfl987tlDermerf2007l:IOstrowskil 
l2002l and references therein). In this paper we, therefore, 
will follow the "jet paradigm", i.e. assume that UHECR par- 
ticle acceleration takes place in jets or is associated with 
them. However, not all AGN have jets, and for those that do 
the jet properties may differ substantially between different 
classes/types (see a discussion in 52). 

There are also quite a few "known unknowns" which affect 
the distribution of the arrival directions of UHECRs and make 
the association with particular classes of objects a non-trivial 
task. Among these are the source distribution, the structure 
of extragalactic and Galactic magnetic fields, the nature of 
sources (transient vs. steady), and the energy spectrum and 
chemical composition at injection (we discuss these topics in 

The Auger Collaboration has reported significant evi- 
dence for anisotropy in the a rrival directions of UHECRs 
(Abra ham et al.ll2007bl l2008ah . The anisotropy signal sug- 
gests a correlation of the e vents with AGN listed in the 
IVeron-Cettv & Veronl J2006) catalog with distances less than 
^100 Mpc, though other sources with a similar distribution 
are not ruled out 3 . The maximum correlation has been found 
for AGN with redshift z < 0.017 (corresponding distance 
D < 71 Mpc), angular separation 9 < 3.2°, and events with 
energy above ~57 EeV. The list of events with energy in ex- 
cess of 57 EeV consists of 27 events, and 20 of these correlate 
with the AGN from the lVeron-Cettv & Veronl d2006l) catalog. 

In this paper we discuss the possible association between 
UHECRs and AGN based on a detailed analysis of a se- 
lected sample of nearby active galaxies contained within the 
9 < 3.2° search radi us of the UHECR eve nts detected by 
Auger, drawn from the lVeron-Cettv & Veronl d2006l) catalog, 
with additional AGN taken from the NASA/IPAC Extragalac- 
tic Database (NED). Our study of the properties of the se- 
lected sample of AGN, described here in detail, shows that 
it consists predominantly of low-luminosity sources such as 
Seyfert galaxies and LINERs, and a handful of radio galaxies. 

In Sj2]we give a brief introduction (not intended as a com- 
prehensive review) into the AGN phenomenology targeting 
non-expert readers; those familiar with the AGN background 
material could safely skip it. In Sj3]we discuss a sample of 
local AGN selected using exactly the same criteria (9 < 3.2° 
and z < 0.018) for each UHECR event that correspond to 
the maximum significance of the correlation reported by the 
Auger Collaboration. We extend the search out to redshifts 
z < 0.037 which doubles the distance to 150 Mpc, and ef- 
fectiv ely covers the range o f horizons for super-GZK parti- 
cles (lAbraham et al.ll2008at lHarari et al.ll2006l) . With such a 
sample, we investigate the properties of the listed 54 active 
galaxies. In particular, we address the two following issues: 
(i) Are the selected objects different in any respect from the 
other numerous local AGNs? (ii) Are they in general able 
to accelerate cosmic rays to ultra high energies? Based on 
the selection criteria, we identify four powerful radio galax- 
ies which could be potential sources of, at least, 7 out of 27 
reported UHECR events. We complete with the discussion 
of the propagation of UHECRs from the sources to the ob- 
server (Sj4]i, and argue that the correlation radius between the 
arrival direction of UHECRs and the sources could be signifi- 
cantly larger than 9 < 3.2° found by the Auger Collaboration 

3 Note that the Auger Collaboration acknowledges the incompleteness of 
the catalog and does not make any serious claims concerning source classes. 



and depends on the UHECR source distribution and the as- 
sumed model of intergalactic and Galactic magnetic fields. In 
Sj5]we summarize our reasoning that the observed correlation 
is most probably chance coincidence given the high density 
of low power AGN in the local universe. 

2. PHENOMENOLOGY OF LOCAL AGN 

A number of weak AGN in the local Universe can be ob- 
served and even resolved due to their proximity which allows 
a broad range of AGN activity to be studied, while only bright 
ones can be seen at large distances. At the current state of our 
knowledge, however, not all AGN types could be recognized 
as potential sources of UHECR. A brief introduction into the 
types of AGN activity is, therefore, warranted as a basis for 
further discussion. 

2.1. Nuclei 

Supermassive black holes (A^bh ~ 10 6 Mq-10 9 Mq) are 
found in a number of galaxies, and there is growing consen- 
sus that all galaxies contain accreting black holes at their cen- 
ters. If the emission of the accreting/circumnuclear matter is 
pronounced, the galaxy is classified as active. Such emission 
includes a non-stellar blue continuum due to the accretion 
disk and strong, but narrow, forbidden emission lines resulting 
from photoionization of the surrounding gas by the disk radi- 
ation. It has been shown that about 85% of all galaxies (and, 
in particular, almost all nearby late-type galaxies) possess de- 
tectable emissio n-line nu clei, which therefore could be clas- 
sified as AGN (Ho 2008). In most cases, however, the nuclear 
luminosity is very weak relative to the stellar radiative output 
of the host galaxy 4 . These nearby low-luminosity AGN can 
be divided into Seyfert galaxies (>10% of local AGN) and 
LINERs (^60%), depending on the ionization level and in- 
tensity of the emission lines . For more details we refer the 
reader to the recent review by Ho (20081). 

The rest of the local AGN assemblage is populated 
by r adio galaxies (mostly low-power ones of the FR I 
type: iFanaroff & Rilevlll974 and BL Lacertae objects (BL 
Lacs). In these sources, the radiative output of the accret- 
ing/circumnuclear gas (which may, spectroscopically, resem- 
ble closely the Seyfert or LINER types) is dominated, or 
at least substantially modified, by a broad-band non-thermal 
emission produced by relativistic jets emanating from the ac- 
tive centers. Also, radio galaxies and BL Lacs are hosted 
almost exclusively by giant elliptical and SO galaxies, while 
the majority of Seyferts and LIN ERs (^90%) are associated 
with s piral galaxies (e.g ., iBahcall et alJ 1199 7: Malka net al.1 
1998; Mar tel et al.lll999l) . Therefore, Seyferts/LINERs and 
FR Is/BL Lacs are morphologically different and this may 
have important consequences/implications for their capabil- 
ity of accelerating particles to the highest observed energies, 
as discussed in the next section. 

Some fraction of AGN exhibit in addition broad permit- 
ted emission lines (FWHM ~ 10 3 - 10 4 kms _1 ) in their 
spectra; the presence/absence of such lines divides further the 
AGN population into type 1/type 2 classes. In the framework 
of the "AGN unification scheme", these two classes are in- 
trinsically identical, and differ only in the orientation of the 

4 If the nuclear emission outshines the starlight of the host galaxy, the 
source is classified as a quasar. 

5 We note that the starburst nuclei of many galaxies, as well as their H II 
regions, also possess ionization lines, which are, however, weaker than the 
ones observed in Seyferts or LINERs, and are photoionized exclusively by 
the starburst activity. 
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accre tion disk symmetry axis to the line of sight dAntonuccil 
119931) . Namely, type 1 AGN are thought to be predomi- 
nantly observed at small inclination angles (roughly <40°), 
while type 2 AGN are viewed at larger angles through a high- 
column-density gas and dust concentrated at pc-scale dis- 
tances from the active center in a torus-like structure (co- 
planar with the inner accretion disks). As a result, the disk 
continuum and broad emission lines produced very close to 
supermassive black holes are heavily obscured in type 2 AGN, 
but not in type 1 sources. The narrow line emission, originat- 
ing at larger distances (> pc) from the active center, is not 
subjected to strong obscuration. 

Such a unifying picture is supported by several observa- 
tional findings, including detections of strong mid-infrared 
(MIR) emission due to the dusty obscuring tori in basically 
all different classes and types of AGN. In addition to the 
MIR/optical/UV non-stellar nuclear radiation, Seyferts and 
LINERs show characteristic X-ray emission extending from 
0.1 keV up to 100 keV photon energies with photon indices 
Tx < 2.0, which is probably prod uced within the hot coro- 
naeof the accretion disks (see, e. g JSvenssonl 1 996t IZdziarskH 
1999). The low-energy (< 2 keV) segment of this continuum 
is typically obscured in type 2 sources, which agrees with the 
unification scheme. We also note that in the framework of the 
unification paradigm, all BL Lacs are considered to be simply 
low-luminosity radio galaxies (FR Is) viewed at very small 
viewing angles (< 10°) to the jet axis. 

2.2. Jets 

Acceleration of particles to ultrarelativistic energies in 
AGN is observationally confirmed (through detection of the 
intense non-thermal emission) to be associated with colli- 
mated fast jets, which are produced by the active black 
hole/accretion disk systems 6 , or in shocks created by a prop- 
agating jet. The strong magnetization, extremely low den- 
sity, non-stationary relativistic and supersonic flow pattern, 
and fin ally turbulent characte r of such jets (see, e.g., the re- 
view by Begelm an et alJl984h . is expected to result in the for- 
mation of non-thermal particle energy distributions extending 
up to the highest accessible energies. Indeed, the very high 
energy 7-ray emission detected from blazars 7 , demonstrates 
the ability of the nuclear (sub-pc scale) relativistic AGN out- 
flows to effici ently accelerate electrons up to 1 — 100 TeV 
energies (e.g JCelotti & Ghise llini 2008). In the case of large 
(kpc)-scale jets, as observed in powerful radio galaxies and 
radio-loud quasars, it has been speculated that proton ener- 
gies in the 10 — 100 EeV range can be re ached as well (e.g., 
iLyufikov & Ouyedl2007tlbstrowskilll998l) . 

Not all AGN are jetted. Also, it seems that only some frac- 
tion of nuclear jets can remain relativistic and well-collimated 
during the propagation through the dense environments of the 
central parts of host galaxies. These issues are particularly 
relevant for Seyferts and LINERs, which differ substantially 
from the bona fide jet sources like radio galaxies and radio- 
loud quasars, lacking in particular well-defined large-scale jet 
structures. We also note in this context, that no galaxy clas- 
sified as Seyfert or LINER has been detected so far at 7-ray 

6 Note, that there is no observational evidence for a significant popula- 
tion of non-thermal particles in AGN accretion disks and disk coronae. In 
fact, hard X-ray/soft 7-ray emission observed from Seyferts is best modeled 
assuming a thermal population of electrons, albeit characterized by high tem- 
peratures of the order of 100 keV. 

7 Blazar class includes BL Lacertae objects and radio loud quasars with 
flat spectrum radio cores. 



(>MeV) energies 8 . 

The jet production and acceleration ef ficiency may depend 
on the morphology of the host galaxy. ISikora et al.l (2007) 
proposed that the difference in the efficiency of the jet pro- 
duction between spiral-hosted AGN (Seyfert galaxies and 
LINERs) and elliptical-hosted ones (low-power FR I radio- 
galaxies, their more powerful analogs classified as FR II ra- 
dio sources, BL Lacs, quasars) may be explained if the spin 
of the central supermassive black holes, J, in spiral galaxies 
is on average lower than in elliptical ones. According to the 
so-called "spin paradigm", the efficiency for the production of 
relativistic jets i n AGN depends on the spin of the supermas- 
sive black hole (Blandford 1990). This has been investigated 
in detail bv lVolonteri et al.l (120071) . who, by means of numer- 
ical simulation, confirmed that indeed the cosmological evo- 
lution of supermassive black holes is expected to lead to such 
a host morphology-related spin dichotomy. If true, this may 
have important consequence for the acceleration of CRs to the 
highest observed energies. Namely, the electromotive force 
of the black hole embedded in an external magnetic field B 
(supp orted by the ac creting matter) is AV ~ Br g (J/J max ) 
(e.g., lPhinnevl[l983l) . where J max = cr g MBK is the max- 
imum value of the black hole spin, and r g is the gravita- 
tional radius of the black hole. Assuming further that the 
magnetic energy density close to the black hole event hori- 
zon is equal to the energy density of the matter accreting at 
the Eddington rate, the maximum energy a test particle with 
an electric charge Ze can reach in such a potential drop is 
£ max ~ZeAl/~3x 1O 2O Z(M BH /1O 8 M ) 1 /2(J/J max ) 
eV. If this picture is correct then many (most) of the lo- 
cal low-power AGN with .Mbh < 10 8 A^q and accreting 
at sub-Eddington rates, and are additionally characterized 
by J <C Jmax as expected in the case of the spiral-hosted 
Seyferts/LINERs, may not have enough potential to acceler- 
ate CRs to ultra high energies jHopkins & Hernquistl l2Q06fc 
ISikora et alj2007llVolonteri et al.ll2007l) . 

Jet activity always manifests itself as non-thermal syn- 
chrotron radio emission produced by ultrarelativistic jet elec- 
trons. Hence, investigation of jet properties are usually ad- 
dressed by detailed radio studies. In the case of local AGN, 
like Seyferts and LINERs, such studies are hampered due to 
the low luminosities and small sizes of the radio structures. 
Additionally, most of these sources, and in principle all lo- 
cal spiral-hosted AGN, show intense star formation activity 
(especially pronounced in the far-infrared [FIR]), see, e.g., 
Ho (2008). Such activity is known for producing the radio- 
emitting outflows due to starburst-driven superwinds which, 
although completely different in origin, may in some (mor- 
phological and spectral) respects closely resemble jet-related 
activity. Thus, care must be taken when making any state- 
ment regarding jet properties in these objects with limited ra- 
dio data. 

The synchrotron/jet origin of the compact cores observed 
in local low-luminosity AGN is support ed by the direct de- 
tections of nuclear jets in these systems (|Ku kula et al1 l!999b 
Middelber g et al.l 120041: lUlvestad etail 1 1 9981) . Interestingly, 
observations of the proper motions of such structures in- 
dicate sub-relativistic bulk velocities on sub-pc/pc-scales, 

8 This may change with the next generation of 7-ray satellites 
(Fermi/Large Area Telescope [Fermi/LAF]) and ground-based Cherenkov 
Telescopes (such as the next phases of HESS, Major Atmospheric Gamma- 
ray Imaging Cherenkov telescope [MAGIC], or Very Energetic Radiation 
Imaging Telescope Array System [VERITAS]). 



4 



Moskalenko et al. 



v < 0.25c (Mi ddelberg et alj|2004t lUlvestad et al.lll999albl) . 
Note that the observed brightness temperatures of the ra- 
dio cores in Seyferts and LINERs, as well as their moderate 
year-timescale variab ility, do not require relativistic beaming 
dMundell et al.l 12000). Also, the observed one-sidedness of 
nuclear jets in these systems is best explained as resulting 
from free-free absorption of radio photons on the surround- 
ing gaseous disks, and n ot due to relativistic (Doppler) effects 
( Middelberg et al. 2004). This constitutes a clear difference 
with the established properties of radio galaxies and quasars. 

We also emphasize that the distribution of nuclear 
jets in Seyferts and LINERs is ra ndom with r espect 
to the host g a laxy stellar disks (iKinney et al.l 120001: 
Nagar & Wilsonl 11999b iPringle et all 119991; ISchmitt et alj 
200 U 120021: iThean et al.ll2001l) . Moreover, misalignments be- 



phology (ellipticals vs. spirals) was found dMoreanti et al 



1999; 



NagaretaD[T999; Rush et al.ll 19961: lUlvestad & W ilso: 



on 



1989). However, more recently it has been suggested that 



tween sub-pc/pc-scale jets and kpc-scale radio structures are 
common in Seyferts, bei ng much larger than those observed 
in other radio-loud AGN dCorbert et alj 1995 IGallimore et alj 
120061: iKharb et alJl2006t llVliddelberg et alJl2004l) . The mis- 
alig nment angle distribution is flat over the whole range 0° to 
90° dMiddelber g et alj2004l) . This may suggest that kpc-scale 
radio structures are not powered by the jets, but rather origi- 
nate in starforming regions. Such an interpretation is some- 
times supported by the fact that the extended radio-emitting 
halos in some Seyferts and LINERS are aligned with the 
galaxy disks, and that their radio powers correlate with the 
FIR fluxes as observe d in regular spiral and starforming galax- 
ies (IThean et alJIfOOlh . 

However, in most of the cases (> 45%) the observed kpc- 
scale radio structures of Seyferts and LINERs do not mor- 
phologically match those of galaxy disks or starforming re- 
gions, show ing also excess radio emission over the ra dio-FIR 
correlation (iGallimore et al.l 120061: IKharb et al.ll2006h . Thus, 
such structures are believed to be powered directly by jets, 
which are, however, substantially different from the ones ob- 
served in other radio-loud AGN. In particular, there is grow- 
ing consensus that Seyferts and LINERs are characterized by 
short-lived low-power jet activity epochs, with jet lifetimes 
tj < 10 5 yrs and jet kinetic luminosities Lj ~ 10 41 — 10 43 erg 
s~ 1 , possibly triggered by minor accretion episodes, which re- 
peat every ~10 6 yrs in different (random ) directions oyer the 
whole Seyfert-type activity epoch < Gyr dCapetfi et alJ[T"999t 
IGallimore et allioOdlKharb et aljfeoOdlSandersI 19841) . 

2.3. Radio activity 

Radio surveys of local Seyferts and LINERs provide fur- 
ther insight into the nature of their activity and reveal typi- 
cally complex, multi-component radio structures, consisting 
of compact unresolved or slightly resolved cores (detected at 
centimeter wavelengths in most of these objects), linear jet- 
like features (observed in ~30% of the local low-luminosity 
AGN), and spherical or elongated diffuse halos/radio bubbles 
(present in, again, ~30% - 50% of the Seyfert and LINER 
population; e.g.. iBaum et al.l [l99l IColbert et alj 119961: 
Galli more et all 120061: iHo & Ulvestadl 120011: iMorganti et alJ 
19991: iRush et alj|1996l: lUlvestad & Wilsonlll989l) . The to- 
tal radio power of these objects at centimeter wavelengths 
is in the range ~ 10 35 — 1 42 erg s" 1 (with me dian value 
£5 GHz ~ 2 x 10 37 erg s'^ lHo & Uivestadll2001l) . which is 
typically slightly higher than the total radio luminosities of 
"regular" spiral galaxies. 

For the compact cores of Seyfert and LINERs, previ- 
ously no dependence in radio power or spectral properties 
on the AGN type (type 1 vs. type 2) or host galaxy mor- 



the small fraction of those Seyferts which are elliptical- 
hosted, and possibly also type 1 Seyferts, a re characterized by 
relatively stronger nuclear radio activity ([Nagar et al.] 120051 : 
IThean et alJl200lt lUlvestad & Holl200ll) . Nevertheless, it is 
established that flat-spectrum radio cores can be present in all 
types of local low-luminosity AGN, and therefore (unlike in 
the case of radio galaxies and quasars) cannot be used as a 
good proxy for the source inclination. 

It has been speculated that the flat-spectrum radio cores of 
Seyferts may not be due to jet activity, but rather some other 
processes like a nuclear starburst, or the em ission of the ac- 
cretio n disks/obscuring tori themselves (e.g.. IGallimore et"aT] 
119971) . However, the relatively high brightness temperatures 
of the Seyfert radio cores, 10 7 K < T b < 10 9 K, points 
to a synchrotro n (and therefore jet-rel a ted) origin of t he ra- 
dio emission dMiddelberg et alJ bool iMundell et all 120001: 
Nagar et alJ 120051) . This is supported by the excess radio 
emission of the local low-luminosity AGN with detected ra- 
dio cores over the radio-FIR correlatio n established for reg- 
ular spirals and starforming galaxies dHo & Ulvestadl 1200 lb 
iRov et al] |T998) which is expected to hold if the radio emis- 
sion i^_ex£hisjy^ly/pred^ activ- 
ity ((CondonitiiJlSol EiiouetiDIilll). Also, recently 
found scaling relations between the nuclear radio emission 
of Seyferts and LINERs and their accretion power (approx- 
imated by either the 2 — 10 keV luminosity of the disk coro- 
nae, non-stellar optical magnitudes of the active nuclei, or Hp 
or O III line luminosities), indicate a strong link between the 
radio production efficiency and the accretion disk parameters, 
as expected in the case of a jet origin for the radio em ission 
dHol2002t IHo & Pengll200ll iPanessa et al]l2TMl2007h . 

3. AGN SAMPLE 

We collected AGN which are possibly associated with the 
UHECR events by search ing both the NASA Extrag alactic 
Database (NED) 9 and the IVeron-Cetty & Veronl d2006l) cat- 
alog. We also extend the search out to redshifts z < 0.037 
(<150 Mpc) which doubles the distance at which maximum 
significance was found by the Auger Collaboration (z < 
0.018, <75 Mpc); this eff ectively covers the range of horizons 
for su per-GZK particles (Abraha m et al.l [2008at lHarari et~aT] 
2006). 

Table Q] lists the AGN found within the 3.2° search 
radius around each UHECR event reported by Auger 
(Abra ham et al.l2008al) . divided into two groups: those within 
redshift range z < 0.018 and 0.018 < z < 0.037 (corre- 
sponding approximately to the luminosity distances gJl < 75 
Mpc and 75 Mpc < d L < 150 Mpc, respectively, Table 
In total, we have selected 54 active galaxies, 27 per each red- 
shift bin. Note that six UHECR events lack any selected AGN 
counterpart located within 75 Mpc, and for two of them we 
also did not find any possible AGN association up to a dis- 
tance 150 Mpc. Four of these are located at low Galactic 
latitudes, \b\ < 12°. Obviously, most of the events possess 
multiple AGN "counterparts" within the assumed search ra- 
dius and redshift range. 

We emphasize that the IVeron-Cetty & Veronl d2006l) and 
NED assemblages are not complete AGN catalogs. Addi- 
tionally, we have found several inconsistencies between these 

9 http : / / nedwww . ipac . caltech . edu/ 
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UHECR Event AGN within 3.2° search radius 

Event number Galactic Coordinates < z < 0.018 0.018 < z < 0.037 







Name 


<?['] 


Name 


e ['] 


25 


(-21.8,54.1) 


NGC 5506 


38 


UM653 


136 










UM654 


138 










UM625 


171 


18 


(-57.2,41.8) 






ESO 575-IG016 


55 










MCG-03-32-017 


191 


26 


(-65.1, 34.5) 






TOLOLO 00020 


133 


2 


(-50.8, 27.6) 


NGC 5140 


116 


2MASX J13230241-3452464 


55 










TOLOLO 00081 


95 










ESO 444-G018 


144 


21 


(-109.4, 23.8) 


NGC 2907 


94 










NGC 2989 


182 






20 


(-51.4, 19.2) 


NGC 5128 (CenA)* 


54 






17 


(-51.2, 17.2) 


NGC 5128 (CenA)* 


139 










NGC 5064* 


162 










NGC 5244 


164 






8 


(-52.8, 14.1) 


NGC 5064* 


47 










IRAS 13028-4909 


118 










NGC 4945 


121 






5 


(-34.4, 13.0) 


IC4518A 


66 






1 


(15.4, 8.4) 






1RXS J174155.3-121157 


106 










HB91 1739-126 


122 


14 


(-52.3, 7.3) 






WKK 2031 


83 


23 


(-41.7, 5.9) 


WKK 4374 


167 






3 


(-49.6, 1.7) 


DZOA 4653-11 


40 










PKS 1343-60 (Cen B) 


41 






27 


(-125.2, -7.7) 










11 


(-103.7, -10.3) 










13 


(-27.6, -16.5) 


ESO 139-G012* 


109 










AM 1754-634 NED03* 


191 






4 


(-27.7, -17.0) 


ESO 139-G012* 


139 










AM 1754-634 NED03* 


166 






10 


(48.8, -28.7) 


CGCG 374-029 


183 


PC 2055+0126 


125 










Mrk 510 


152 


19 


(63.5, -40.2) 


PC 2207+0122* 


98 


PKS 2201+044 


109 










NGC 7189* 


151 


7 


(58.8, -42.4) 


PC 2207+0122* 


189 


NGC 7189* 


108 


16 


(-170.6, -45.7) 


NGC 1358* 


51 


Mrk 612 


81 






NGC 1346* 


66 


Mrk 609* 


133 






NGC 1320 


130 


KUG 0322-063A* 


137 


12 


(-165.9, -46.9) 


NGC 1346* 


152 


Mrk 609* 


169 






NGC 1358* 


166 


KUG 0322-063A* 


171 


15 


(88.8, -47.1) 


NGC 7626 


86 


2MASX J23274259+0845298 


122 






NGC 7591 


186 


NGC 7674 


125 


24 


(12.1, -49.0) 


NGC 7130 


117 


6dF 32 132022-334254 


127 






NGC 7135 


126 






22 


(-163.8, -54.4) 


NGC 1204 


96 


MCG-02-08-039 


174 


9 


(4.2, -54.9) 


IC5169 


131 


ESO 404-IG042 


76 










PKS 2158-380 


126 


6 


(-75.6, -78.6) 


NGC 0424 


25 


ESO351-G025 


153 










ESO 352-G048 


163 



NOTE. — 6 denotes the separation between an UHECR event and an AGN. Stars denote A GN falling within the 3.2° search radius of two different events. 
AGN marked with bold font are those which do not appear in the Veron-Cetty & Veron 1 2006) catalog. 



two databases regarding source classifications. This results 
in a few sources showing some (weak) level of AGN activ- 
ity (according to one but not the other catalog) which might 
fulfill the criteria for a possible association with the detected 
UHECR events, but these are not included in our dataset 10 . 
However, a few other objects which do not obviously possess 
an active nucleus (but rather only H II central activity, like 
NGC 4945, NGC 5244, NGC 2989, NGC 7135, and espe- 

10 For example, we did not include the object ESO 383-G18 in our dataset, 
since its more accurate position from 2MASS as given in NED makes it 
192.2' a way from the location of CR #2. With the less accurate position 
listed by Veron-Cetty & Veron 12006), however, this source is <192' from 
the considered CR event. 



cially nearby IRAS 13028-49 1 ') are included and classified as 
"possibly LINERs." Furthermore, the regions around 3 events 
(#7, #1 2, #16) are cover ed in the Sloan Digital Sky Survey 
(SPSS: lYork et al] |2000) leading to a somewhat misleading 
higher density of known AGN in these regions in compari- 
son to other event re gions. We have therefore om itted the 12 
SDSS AGN in the IVeron-Cettv & Vero n (2006) list falling 
within the search volume from Tables Q] [2] and listed them 
separately in Table[3] These objects, all classified as Seyferts, 
are all located at redshifts z > 0.018 and display similar bulk 
properties to the Seyferts we have included in our discussion. 
Most of the selected objects (50 out of 54) are relatively weak 

See IStrauss et all (1993) . 
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Table 2. Properties of selected AGN 



Name 


Type 


RA 


DEC 


z 




UHECR 


e 


Sr 


Ref. 




Sx 


Ref. 






(J2000.0) 


(J2000.0) 




[Mpc] 




['] 


[mJy] 




[Jy] 


[cgs] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 5506 


Sy2 


1 A L. 1 1 1 A O - 

14nl3ml4.8s 


at j 1 n„ 

— 0idlzmz7s 


0.0062 


29 


(—21.8, 54.1) 


38 


227 


G06 


8.4 


12 


S06 


T T1V /I CZ1 

UM 6j3 


Sy2 


14nloml j.js 


— Oldz /mjjs 


0.0365 


IDO 


(—21.8, 54.1) 


koo 


0.9 C 


B95 


n 1 
U.4 






T T1V /I CZA 

UM 6j4 


Syz 


1 A U 1 Am 1 O In 

14nl6miy. /s 


— Oldzjmlos 


0.0369 


160 


(—21.8, 54.1) 


138 


^-n oc 
<U.2 


B9J 








UM 625 


Sy2 


14n00m40.6s 


A 1 JCC 1 Cl 

— Uld55ml5S 


0.0250 


109 


/ TIO C A 1\ 

(—21.8, 54.1) 


171 


l.l c 


B95 


0.3 






hSO 5 /5-IG016 


S2 


12n52m36.2s 


— zla54m4os 


A ATI A 

0.0zz9 


100 


/ C"7 T A 1 0\ 

(—5 l.l, 41.8) 


55 


1.7 C 


C98 








MCG-03-32-017 


L(?) 


12n38m00.5s 


— z0a07m51s 


A AT OA 

U.UzoU 


123 


(—5 /.2, 41.8) 


191 


<0.8 C 


C98 








1 ULULU 20 


Syl 


IznlzmzU.Us 


'fQAA Om /I /if 

— zoa4om4os 


A AT A A 
0.030U 


1 3 1 


(—65.1, 34.5) 


133 


<0.5 C 


C98 








zMAoA J 1323 


L(?) 


I 3n23m02.4s 


1 A A dm A If 

— j4cDzm4 /s 


0.0261 


1 1 3 


\— jU.o, 2/. 6) 


55 


10. 8 C 


C98 








1 OLOLO 8 1 


Sy2 


13nl9m38.6s 


— 33a22m54s 


0.0291 


127 


(—50.8, 27.6) 


95 


<0.8 C 


C98 








NGC 5140 


L(?) 


13n26m21./s 


— jJa5zm0os 


0.0129 


58 


/ cao n /;\ 

( — 50.8, 27.6) 


1 16 


29.2 


S89 


0.8 






£,oW 444-LrOlo 


t /Tt 
L(.') 


1 jnzzmjo.os 


— jza4jin4zs 


n moo 
o.ozyz 


1 T7 
1 2 / 


( en c 07 &\ 
\ — jU.o, Z/.O) 


144 


"7 CC 
/. J 


cyo 








wrjr 7Q07 


t n\ 
H '.) 


OQlyl 1 m3A 7c 
\}y\\ D 1 111 JO. IS 


— 1 UU' 1 t i +lllU J a 


0070 


33 


( 1 OQ d 93 8"! 


OA 


4 4 C 


M92 


4 






NGC 2989 




09h45m25.2s 


— 18d22m26s 


0.0138 


62 


(—109.4, 23.8) 


182 


7.6 C 


C98 


1.7 






NGC 5128 (Cen A) 


FR I 


13h25m27.6s 


— 43d01m09s 


0.0018 


3.4 


(—51.4, 19.2) 


54 


53792 


G94 


162 


30 


B06 














(-51.2, 17.2) 


139 












NGC 5064 


L 


13hl8m59.9s 


-47d54m31s 


0.0099 


45 


(-51.2, 17.2) 


162 


4.0 C 


M03 


3.3 


















(-52.8, 14.1) 


47 












NGC 5244 


L(?) 


13h38m41.7s 


-45d51m21s 


0.0085 


38 


(-51.2, 17.2) 


164 


5.7 C 


M03 


1.9 


— 


— 


IRAS 13028-49 


Sy(?) 


13h05m45.5s 


-49d25m22s 


0.0012 


8.3 


(-52.8, 14.1) 


118 


<22 c 


B99 


6.1 


— 


— 


NGC 4945 


Sy(?) 


13h05m27.5s 


-49d28m06s 


0.0019 


3.8 


(-52.8, 14.1) 


121 


2953 


G94 


359 


0.5 


L04 


IC 4518A 


S2 


14h57m41.2s 


-43d07m56s 


0.0163 


70 


(-34.4, 13.0) 


66 


64.9 C 


M03 





1.9* 


B07 


1RXS J174155 


Syl 


17h41m55.3s 


-12dllm57s 


0.0370 


156 


(15.4, 8.4) 


106 


1.4 C 


C98 




2.9^ 


B07 


HB91 1739-126 


Svl 


17h41m48.7s 


-12d41m01s 


0.0370 


156 


(15.4, 8.4) 


122 


<0.5 C 


C98 








WKK 2031 


S2 


13hl5m06.3s 


-55d09m23s 


0.0308 


133 


(-52.3,7.3) 


83 


154 


G94 


41.0 






WKK 4374 


Sy2 


14h51m33.1s 


— 55d40m38s 


0.0180 


77 


(—41.7, 5.9) 


167 


6.2 C 


B99 




2.0^ 


B07 


\J/,\ J r\ fOJjt-l 1 


ayi 


1 Jll'-t 1 111 J 0. US 


A0H37m0zlc 
— OUUj /lllUT-a 


01 9Q 


SA 

JO 


( AQ A 1 7^ 


-+U 


^3 s c 

J.o 


RQQ 

D77 




7 8^ 


B07 


rrvo 1j4j-6U (Cen d) 


rK 1 


1 11-1/1 tm/IO fir- 
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o.oizy 


J 6 


\~ 4y.o, l.l) 


A 1 

4 1 


77 1 nn 
2 / 1U0 


wyu 




n 7^ 
0.2 


»inc 
MID 
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»yz 


1 7h37m3Q 1 c. 
1 / 11 j / 111 J y . 1 a 


SQHSAm97c 
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7 1 
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1 U 7 
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M03 


7 
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^0 Q c 
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(88.8, —47.1) 


86 
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0.03 


T05 
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zMASX J232/ 


Sy2 
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z3nz /m4z.6s 


1 f\0 A A C O A„ 

+0od45m30s 


A Al A A 

0.0294 


1 1 8 


/O OO A 1 1 \ 

(88.8, —47.1) 


122 


<0.5 C 


C98 








M/T "7/^-7,1 

NGC 76/4 


Sy2 


z3nz/m5o./s 


+08d4om45s 


A AT OA 

0.0289 


1 16 


/O O O A 1 1 \ 

(88.8, —47.1) 


125 


90. 6 C 


C02 


5.6 


0.05 


L04 


JNCC /jVI 


L 


zjnlomlo.js 


+Uod3jm0ys 


0.0165 


64 


(Q Q O A *7 1 \ 

(00.0, —4 /.I ) 




71 A C 

21 .4 


C02 


7.2 






JNGC / 13U 


Sy2 


zln4omiy.Ds 


1 A A C7m(Kn 

— 34dD /itiOjs 


0.0161 


64 


/" 1 1 1 ac\ n\ 
(12.1, — 4y.U) 


1 17 


63.3 


M03 


1 6 


A AA/- 

0.006 


Gm06 


NGC 7135 




9 1 h49m4(S 0s 

£4 1 11^ v 111 iU . \J [a 


— 34H59m35s 


0088 


33 


(12.1 —49.0) 


126 


2.3 C 


C98 


0.2 






6dFJ2 132022 


Syl 


21h32m02.2s 


-33d42m54s 


0.0293 


120 


(12.1^ -49.0) 


127 


l.l c 


C98 








NGC 1204 


L 


03h04m39.9s 


-12d20m29s 


0.0143 


57 


(-163.8, -54.4) 


96 


9.4 C 


C90 


7.8 






MCG-02-08-039 


Sy2 


03hOOm30.6s 


-Ild24m57s 


0.0299 


123 


(-163.8, -54.4) 


174 


3.6 C 


C98 


0.5 






ESO 404-IG042 


S2 


22hl3ml7.5s 


-37dOOm58s 


0.0340 


140 


(4.2, -54.9) 


76 


2.2 C 


C98 


0.7 






PKS 2158-380 


FRII 


22h01ml7.1s 


-37d46m24s 


0.0333 


137 


(4.2, -54.9) 


126 


590 


W90 


0.3 






IC5169 


Sy2 


22hl0ml0.0s 


-36d05ml9s 


0.0104 


39 


(4.2, -54.9) 


131 


8.3 C 


M03 


3.4 






NGC 0424 


Syl 


01hllm27.6s 


-38d05m00s 


0.0118 


46 


(-75.6, -78.6) 


25 


9.6 C 


C98 


1.8 


0.1 


U05 


ESO351-G025 


Sy2 


00h58m22.3s 


-36d39m37s 


0.0346 


143 


(-75.6, -78.6) 


153 


2.4 C 


C98 








ESO 352-G048 


Sy2 


01h20m54.7s 


-36dl9m26s 


0.0322 


132 


(-75.6, -78.6) 


163 


<0.5 C 


C98 









RE FERENCES. — (B95 ) fBecker et all 19951 (BO^rBeclanann et alJ200ft (B0 7) |Bird et allzTOt (B99)|Bock et al|1999l (B94)|Brinkmann et all 19941: (C90)lCondon et aUl99Q| ; 
(C98) Condon et al. 1998; (C02) Condon et al. 2002; (G06) Gallimore et al. 2006; (Gm()6) Gonzalez- Martin et al. 2006; (G94) Gre.sory et al. 1994; (G05) Guainazzi et al. 2005; (L04) 
Lutzetal. 2004; (M05) Marshall et al. 2005; (M03) Mauch et al. 2003; (M92) Mollenhoff et al. 1992; (S89) Sadler et al. 1989; (S06) Shinozaki et al. 2006; (T05) Tajeretal. 2005; 
(U05) Ueda et al. 2005: (W92) White & Becker 1992; (W90) Wright & Otrupcek 1990. 



NOTE. — [1] Name of the source. [2] AGN classification (Syl/Sy2: Seyfert galaxy of the type 1/type 2; L: LINER; FR I: radio galaxy of the FR I type; BL: BL Lacertae object). 
[3-4] Equatorial coordiantes (J2000.0). [5] Redshift of the source. [6] Lu minosity distance for the a ssumed cosmology (Hp — 73 km s - 1 Mpc - 1 , Om — 0.27,Oa — 0.73) except 
for NGC 5128 (Cen A) and NGC 4945, where the distances are known i Karachentsev et al. 2007). [7] Galactic coordiantes for the nearby UHECR event. [8] Separation between an 
AGN and a nearby UHECR event. [9] The total 5 GHz flux in mJy units ( c fluxes conv erted from the ones provided at lower frequencies assuming radio spectral index a — 0.7). [10] 
References for column 9. [ll]The total 60 /im flux in Jy units from the IRAS su rvey i Moshir et al. 199(1). [1 2] The total observed 2 — 10 keV flux in X IP -11 erg cm -2 s" 1 units 
( l fluxes provided at lower photon energy range: 0.5 — 7 keV in Marshall et al. 1 2005), or 0. 1 — 2.0 keV in Brinkmann et al. 1 1994) and Taier et al. 1 2005); ^fluxes converted from 
the ones provided at higher photon energy range 40 — 100 keV in Bird et al. ( 2007), assuming X-ray photon index Tx — 2.0). [13] References for column 12. 
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Table 3. SPSS AGN dHAO et al.I2005I) . 



Event Number 


Name 


Type 


z 


e 


(1) 


(2) 


(3) 


(4) 


['] 

(5) 


1 


oijio jzijzjy,u/-uuuyuj.4 


Syi 


U.Uzo 


1 oo 




SDSS J2205 15.43-010733.3 


Sv2 


0.032 


7 


12 


SDSS J033458.00-054853.2 


Syl 


0.018 


123 




SDSS J034545. 16-071526.8 


Sy2 


0.022 


172 




SDSS J033955.97-063228.9 


Sy2 


0.031 


113 




SDSS J033713.31-071718.0 


Sy2 


0.033 


53 




SDSS J032329.63-062944.1 


Sy2 


0.034 


182 




SDSS J034330.25-073507.4 


Sy2 


0.036 


135 


16 


SDSS J033458.00-054853.2 


Syl 


0.018 


97 




SDSS J033955.97-063228.9 


Sy2 


0.031 


179 




SDSS J033713. 31-071718.0 


Sy2 


0.033 


190 




SDSS J032329.63-062944.1 


Sy2 


0.034 


168 



NOTE. — [1] UHECR Event Number. [2] Name of the source (note that 
source positions are implicit in the SDSS names). [3] AGN classification 
(Syl/Sy2: Seyfert galaxy of the type 1/type 2). [4] Redshift of the source. [5] 
Separation between an AGN and a nearby UHECR event. 



Seyfert galaxies and LINERs, while only a small fraction (4) 
are bright and well established jet sources (three radio galax- 
ies and one BL Lac object). 

It is probable that more Seyferts/LINERs are missing (as 
evident in the SDSS covered regions) and/or mis-classified 
in the presented AGN collection. However, we expect much 
fewer radio galaxies/BL Lacs to be omitted. For an order 
of magnitude estimate, we derive an average number den- 
sity of radio galaxies within the redshift range z < 0.037 
of ~ 3 x 10 -6 Mpc~ 3 (consistent with an estimate of the 
space density of UHECR sources given below) assuming the 
total radio luminosity range between 10 38 erg s _1 and > 10 42 
erg s . This was derived using the 151 MHz luminosity 
function of low -power radio sources (including classical FR I, 
and FR II rad io galaxies with weak/absent emission lines) of 
IWillott et alj d2001l) . converting their model C to the mod- 
ern cosmology adop ted in this paper (discussed in detail in 
IStawarz et all l2006h and the total 151 MHz luminosities to 
5 GHz ones assuming a radio spectral index of a# = 0.7. 
The derived number density of local radio galaxies is there- 
fore about three orders of magnitude lower than the number 
density of Seyfert galaxies (see below), resulting in a total of 
—42 expected radio galaxies within the co-moving volume of 
the local Universe V = 0.014 Gpc 3 . We note that with the 
assumed minimum radio luminosity one order of magnitude 
lower, i.e. 10 37 erg s — 1 , the number density of radio galaxies 
increases by a factor of —4. At such low radio l uminosities, 
howev er, the luminosity function constructed by lWillott et al.1 
d2001l) is not well determined. 

With the small number of nearby radio galaxies, we expect 
that the majority of these objects have already been identi- 
fied as such through radio imaging studies. The luminosities 
of the radio galaxies/BL Lac in our sample (down to 4 x 10 39 
erg s _1 , the 5 GHz luminosity of Cen A), correspond to a rela- 
tively bright flux limit (—100 mJy source 12 at D = 160 Mpc). 
The morpholog ies of such bright radio sources hav e been well 
surveyed (e.g. JOwen et al.ll995USadler et al.ll989l and refer- 
ences therein) although admittedly, there is sparser coverage 
in some parts of the sky, in , e.g., the Southern hemisphere 
(Burg ess & Hun stead 2006). The Galactic plane is particu- 

12 For aj = 0.7, this corresponds to ~250 mJy at 1.4 GHz. 



larly susceptible to having unidentified local radio galaxies 
but the radio imagin g surveys toward the plane are maturing 
dHelfand et al.ll20 06) and hard X -ray surveys are b eginning to 
help rectify the situation (e.g., Mol ina et al.ll2007l) . 

Table|2]lists the main properties of all the selected AGN col- 
lected from the literature, in particular the total radio fluxes 
at 5 GHz (or upper limits for these), IRAS fluxes at 60 fim 
(if available), and 2 — lOkeV fluxes (again, if available). 
In some cases, the flux conversion from values provided at 
lower radio frequencies (or in higher X-ray photon energy 
range) was necessary, and we performed this assuming typ- 
ical values of the radio spectral index cur = 0.7 (X-ray pho- 
ton index Fx = 2.0). Note that 12 sources listed in Ta- 
ble |2]have only upper limits for the total radio emission which 
were derived by us fr om the FIRST { Becker et al. 1995), NVSS 
(ICondon et al.ll 19981) . or SUMMS (Boc k et al.lll999h surveys, 
and most of them are indeed very weak radio sources (<mJy 
level). 

Only in the case of the four selected FR I/BL Lac sources 
(Cen A, Cen B, PKS 2158-380, PKS 2201+044) is the pres- 
ence of relativistic jets - as ones considered in different mod- 
els for the acceleration of UHECRs - certain. As for the rest 
of the selected AGN, one should ask if jet activity can be as- 
certained and, if so, what are the properties of the jets in these 
sources. Are these jets relativistic? What are their inclina- 
tions? Do they differ somehow from the jets observed in other 
Seyfert and LINER galaxies? Such questions are relevant, 
because the space density of the local low-luminosity AGN, 
just like the ones selected in our sample, is high. For exam- 
ple, the space density of bright galaxies (—22 mag < Mb < 

— 18 mag) showing S eyfert activity is ~1.25 x 10~ 3 Mpc~ 3 
dUlvestad & Holl2001l) 13 which is -10 - 10 3 times larger than 
the space density of UHECR sources, — 10~ 6 - 10~ 4 Mpc -3 , 
derived from a comparison of simulation of particle propaga- 
tion in the local univers e and the Akeno Giant Air Shower 
Array (AGASA) data dBlasi & De Marcol I200H ISigl et al l 
120041: iTakami et al.ll2006l) . Note that iDubovs kv et al. (l2000h 
gives an estimate of the space density of UHECR sources 
—6 x 10~ 3 Mpc~ 3 based on 7 observed doublets (only two 
events have estimated energy >100 EeV) and assuming the 
distance <25 Mpc for protons >100 EeV. When rescaled to 
the distance <75 Mpc, it gives — 2 x 10~ 4 Mpc~ 3 , con sis- 
tent with other estimates. Recently Takami & Satol d2008h es- 
timated a space density of UHECR sources — 10~ 4 Mpc~ 3 
based on the Auger data. Thus, restricting the investigations 
to redshifts z < 0.037, corresponding to the luminosity dis- 
tance c?l < 156.4 Mpc and co-moving volume V = 0.014 
Gpc 3 , there are —1.8 x 10 low-luminosity AGN with sur- 
face density (if distributed isotropically over the whole sky) 
—1.4 x 10 3 per steradian. In other words, there should be 

— 14 AGN within the search radius 3.2°. The much lower rate 
of our possible identification given in Table[T]is due to incom- 
pleteness/lack of the AGN surveys in different parts of the 
sky, and especially within the Galactic Plane (see a decrease 
in the identification rate for Galactic latitudes \b\ < 12° in 
Table Q]). Therefore, if the selected Seyferts and LINERs do 
not differ from the other local AGN of the same types, the 
claimed correlation between the UHECR events and local ac- 
tive galaxies should be considered as rather unlikely, resulting 
from a chance coincidence, if the production of the highest 

13 This is about ~10 3 times larger than the local space density of bright 
quasars, and ~10 times sm aller than the space de nsity of "regular" galaxies 
with comparable brightness lUlvestad & Ho 200T1). 
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FIG. 1. — (a) Number distributions of 5 GHz total radio powers for 
Se yfert galaxies d e tected at radio frequencies from the sample constructed 
by Ho & Ulvestad 12001, red solid), Seyfert/LINER galaxies from our sam- 
ple detected at radio frequencies (green dashed), and radio galaxies from our 
sample (blue dotted), (b) Number distributions of upper limits to 5 GHz to- 
tal radio powers for S eyfert galaxies d etected at radio frequencies from the 
sample constructed by Ho & Ulvestad 12001, red solid), and Seyfert/LINER 
galaxies from our sample detected at radio frequencies (green dashed). The 
histograms are slightly shifted relative to each other in horizontal direction 
for clarity. 



energy CRs is not episodic in nature, but operates in a single 
object on long (> Myr) timescales. In addition, if the selected 
sources do not show significant jet activity, there are no rea- 
sons for expecting them to accelerate CRs up to 10-100 EeV 
energies. 

A brief summary of the radio properties of local low- 
luminosity AGN given in ij2] allows us to conclude that the 
particular Seyferts and LINERs selected in Table |2] are most 
likely jetted, but do not differ from the other analogous 
sources of the same type. In particular, complex radio mor- 
phologies consisting of compact cores, one-sided jets, and 
extended halos, as found in several of the selected objects, 
are typical for the Seyfert/LINER-type jet activity (see Sj2). 
Such jets are expected to be sub-relativistic (v < 0.25c on pc- 
scales), low -power (Lj < 10 43 erg s _1 ), precessing and short- 
lived (tj < 10 5 yrs). Indeed, 5 GHz powers for the selected 
Seyferts and LINERs, being in a range L 5 GHz ~ 10 37 — 10 42 
erg s _1 , are typical for the other AGN of the same kind, 
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FIG. 2. — X-ray (2 — 10 keV) luminosity versus 5 GHz total luminos- 
ity for Seyfert galaxies from the sample constructed by Panessa et al] 120071 
black circles), Seyfert/LINER galaxies from our sample with the X-ray and 
5 GHz fluxes provided (magenta crosses), radio galaxies from our sample 
with t he X-ray fluxes prov ided (open squares). In the case of the sam- 
ple of |Panessa et al. ( 2007), only sources detected at X-ray and radio fre- 
quencies were considered. Red solid line indicates the best fit log Lx = 
0.97 logLfl + 5.23 for Seyfert galaxies, and blue dotted line denotes the 
best fit \ogLx = 0.97 log L r + 2.42 for low-luminosity radio galaxies, 
both as given by Paness aet all 



and are significantly lower than the radio powers of 4 se- 
lected radio galaxies/BL Lac. The median values of these, 
r^2 x 10 38 erg s _1 (including only radio-detected sources), 
se ems to be higher than the appropriate median values given 
bv lHol d2008l) . but this may be simply due to selection effects. 
In fact, all 12 objects for which only upper limit regarding the 
radio fluxes are provided, have Lsghz < 5 x 10 37 erg s _1 . 

Figure QJ shows the number distributions of 5 GHz to- 
tal radio powers for Seyfert galaxies detected at radio fre- 
quencies from t he op tically selected sample constructed by 
iHo & Ulvestadl d2001l) , the Seyfert/LINER galaxies from our 
sample detected at radio frequencies, and radio galaxies from 
our sample. At first glance, our Seyferts/LINERs appear 
systematically more radio luminous in comparison to the 
IHo & Ulve stad (2001) sample by 1-2 dex. However, we note 
that the upper limits we obtain from the available "all-sky" 
maps (NVSS, S UMMS) are s ystematically higher than the 
ones obtained bv lHo & Ulvestadl from their pointed observa- 
tions (FigureQJ) also by 1-2 orders of magnitude, making it 
difficult to compare the two distributions directly. The radio 
galaxies are characterized by larger radio luminosities than 
Seyferts/LINERs, as expected. 

The range of X-ray luminosities (if available) of objects in 
our sample is in the range L 2 _iokcV ~ 10 40 — 10 44 erg s _1 
(with the median ~ 3 x 10 41 e rg s _1 again slightly higher than 
the one provided by Ho 2008), which is comparable to the typ- 
ical 2-10 keV lumin osities observed in local low -power AGN 
dPanessa et al.l l2~007). Moreover, the logarithm of the ratio of 
the X-ray and radio luminosities, log(L5GHz/-^2— iokev)> al- 
though widely scattered in a range between (—1.3) and (—5.6) 
with median (—3.2), is in agreement with the values foun d in 
other local low-luminosity AGN (see lPanessa et al.ll2007l) . 

Figure [2] shows the X-ray (2 — 10 keV) luminosity versus 
5 GHz total lu minosity for Seyfert galaxies from the sample 
constructed bv lPanessa et alj (|2007), Seyfert/LINER galaxies 
from our sample with the X-ray and 5 GHz fluxes provided, 
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FIG. 3. — 60 /im luminosity versus 5 GHz tota l lumi nosity for Seyfert 
galaxies from the sample constructed by Roy et al. 1 1998, detections: black 
circles, upper limits: open circles), Seyfert/LINER galaxies from our sample 
with the 60 (im and 5 GHz fluxes provided (magenta crosses), radio galaxies 
from our sample with the 60 fim flu xes provided (open squares). In the 
case of the sample of Roy et al. 1 1998), sources classified as quasars were 
omitted, and radio luminosities provided at 2.4 GHz were converted to 5 GHz 
luminosities assuming radio spectral index an = 0.7. 



and radio galaxies from our sample with the X-ray fluxe s pro- 
vided. In the case of the sample of Panessa et al. (2007), only 
sources detected at X-ray and radio frequencies were consid- 
ered. The bulk of our selected Seyferts/LINERs match well 
the Lx — Ln correlation established for nearby Seyferts. Our 
selected radio galaxies and a few Seyferts are over-luminous 
in the radio for a given X-ray luminosity when compared to 
the Seyferts. However, they match well the Lx — Lr corre- 
lation esta blished for nearby "low-luminosity radio galaxies" 
(LLRGs. lPanessa et al.ll2007l) . 

Finally, the ratio of FIR and radio luminosities for the se- 
lected Seyferts and LINERs agrees with what is observed in 
other analogous sources. Namely, the median value in our 
sample is log(i60^m/^5GHz) ~ 5.37 for Seyferts and LIN- 
ERs. These ca n be c ompared with the medians claimed by 
iHo & Ulvestadl (1200 ll) for Seyferts, -5.31, and for regular 
spiral galaxies, —5.64, implying that the Seyferts and LIN- 
ERs included in our sample are not more than twice brighter 
in radio than expected if all the radio emission is due to star- 
forming activity, in agreement with the established properties 
of other Seyferts. We note that the analogous ratios for the 
two radio galaxies included in the sample and detected at FIR 
are significantly lower, log(L6o M m/£5GHz) ~ 2.8 and 3.5, as 
expected. 

Figure [3] shows the 60 /im luminosity versus 5 GHz total 
lu minosity for Seyf ert galaxies from the sample constructed 
by iRov et alj {1998), Seyfert/LINER galaxies from our sam- 
ple with the 60 fim and 5 GHz fluxes provided, and radio 
galaxies from our sample w ith the 60 urn fl uxes provided. 
In the case of the sample of IRov et al.l (1 19981) . sources clas- 
sified as quasars were omitted, and radio luminosities pro- 
vided at 2.4 GHz were converted to 5 GHz luminosities as- 
suming a radio spectral index an = 0.7. All our selected 
Seyferts/LINERs match the FIR-radio correlation established 
for other nearby Seyferts. Our selected radio galaxies are 
over-luminous in radio for a given FIR luminosity with re- 
spect to Seyferts, as expected. 



3.1. Selected radio galaxies 

As already discussed in $2] radio activity in low-power 
active galaxies of the Seyfert or LINER type differ sub- 
stantially from that observed in well established jet sources 
like radio galaxies and radio-loud quasars. Indeed, there 
are some important reasons for such a difference. In par- 
ticular, Seyferts and LINERs are usually hosted by spiral 
(disk) galaxies, while radio galaxies and radio-loud quasars 
are typically hosted by giant ellipticals. It was noted re- 
cently that merger episodes triggering jet activity (by shap- 
ing the accretion processes and even determining spins of 
supermassive black holes), proceed differently depending on 
the host properties dHopkins & Her nquist 2006; Siko ra et ail 
2007; IVolonteri et aU\200% ~ 

Of the 54 selected AGN (Table |2), 3 are classified as radio 
galaxies and one is a BL Lac object, i.e., under the unifica- 
tion sc heme, it is a FR I obse rved with a small jet viewing 
angle ( lUrrv & Padovani| [l995). This group includes Cen A 
(Figure [4), which is characterized by a well-known FR I ra- 
dio morphology with a one-sided (—4 kpc lon g) jet and gi- 
ant (-0.5 Mpc) radio lobes (e.g., Ilsraellll998l). This is the 



only source in the sample detected in 7-rays ( Sreekum ar et alj 
119991: ISteinle et alj Il998t> . Radio maps of Cen B reveal a 
one-sided well-defined FR I jet extending from pc- to kpc- 
scales, and an edge-brightened structure on the o pposite lobe 
more characteristic of a powerful FR II (Figure 131 iJones etaTI 
12001 . From our new VLA map of PKS 2158-380 (Fig- 
ure its radio morphology is characteristic of a FR II with 
bright compact features at the outer edges of its lobes. It is, 
however, relatively underluminous for a FR II (Li.4GHz ~ 
4 x 10 24 W Hz -1 ) and can be considered an intermediate 
object. PKS 2201+044 is classified as a BL Lac object with 
an asymmetric core-jet radio st ructure (Figure U}, an d an ex- 
tended (—100 kpc) radio halo (lAugusto et al.lfl998l) . These 
four objects are quite representative of the local population of 
radio galaxies. Note that Cen A is exceptionally bright and 
extended in the sky only due to its proximity. 

In fact, due to its large extent on the sky, Cen A may be con- 
sidered as being likely associated with more than two Auger 
events (see Figure |4). Note that Table Q] indicates a possible 
association of Cen A with only two CR events, because the 
3.2° separation of the UHECR event from the Cen A nucleus 
is assumed when compiling the list of the selected AGN). 
Considering the 9 events plotted in Figure [4] around the giant 
radio structure of Cen A, their detection rate appears steady 
(almost every third event detected by Auger, except for the 
larger gap between events #8 and #14). If the rate is indeed 
steady, as would be expected if giant Mpc-scale radio lobes of 
Cen A are the acceleration sites of these events, this will be 
easily tested with additional 2-4 years of Auger observations. 

Considering the four events closest to Cen A, the positions 
of these four events are roughly aligned with the axis of the ra- 
dio lobes, whic h is also aligned with t he super-galactic plane 
(cf. Figure 2 of lAbraham et al.ll2008af) . The latter two events 
are closest to the center of Cen A and the former two are co- 
incident with other AGN in the field. It is of interest to note 
that one of the field AGN is NGC 4945 (event #8), which 
is located in the Centaurus group (i.e., it is at the same ap- 
proximate distance as Cen A). NGC 4945 hosts a less power- 
ful radio source than Cen A and is dominated by extended 
emission from the g alaxy and a compact non-thermal core 
dElmouttie et al.lll997l) . Also, of the 6 UHECR events without 
an associated nearby (z < 0.018) AGN within a 6 < 3.2° cir- 
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FIG. 4.— Radio map (at 408 MHz from IHaslam et all 19821) of the 35° X 
35° field centered on the nearby radio galaxy Cen A. The total extent of the 
north-south radio lobes is ~9° and is centered on the AGN (the bright white 
region near the center of the field). The r = 3.2° circ les mark the p ositions 
of the UHECR events detected in the field by Auger ( Abraha met al.H2008al) 
The nu mbers correspond to the event number as provided in Abraham et al. 
12008a), and also in our Table [T] Note event #3 corresponds most closely to 
Cen B, a bright spot near the center of the circle, shown with higher resolution 
in Figure|3] 
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FIG. 5. — Radio image of the 0.75° X 0.75° field centered on the radio 
galaxy Cen B. The 43" resolution image was obtained with the MOST at 
843 MHz by McAdanl jl991l) . The location of the closest CR detected by 
Auger (#3; c.f. Figure |4) is indicated by the arrow pointing away from the 
radio nucleus. 



cular area, 2 are in the plotted field (#14, #26). However, #14 
is the one closest to the Galactic Plane, where it is more diffi- 
cult to identify AGN, and #26 is the one furthest from Cen A. 
As discussed in the next Section a larger deflection angle for 
the events close to the Galactic plane is possible, which could 
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FIG. 6. — Radio image of the 200" X 200" field around the radio galaxy 
PKS 2158-380. This VLA 4.9 GHz image at 2" resolution was made using 
a multi-configuration dataset consisting of a 1 hr observation in Dec 1983 
(program AT45) and 10 min observation from Jun 1997 (AK444). The latter 
dataset was obtained through the NRAO VLA Archive Survey (NVAS). The 
location of the closest cosmic ray detected by Auger (#9) is indicated by 
the arrow pointing away from the radio nucleus. (This NVAS image was 
produced as part of the NRAO VLA Archive Survey, (c) AUI/NRAO.) 
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FIG. 7.— Radio image of the 100" X 100" field around the BL Lac object 
PKS 2201+044 showing its one-sided radio jet. The image is a lower resolu- 
tion (1 " beam) version of the VLA 8.5 GHz data published in Samb rana et alj 
120071) . The location of the closest CR detected by Auger (#19) is indicated 
by the arrow pointing away from the radio nucleus. 



mean that Cen B could be associated with more than 1 event. 

4. PROPAGATION AND COMPOSITION OF UHECRS 

The propagation of the UHECRs from the sources to the ob- 
server is not rectilinear due to deflection by intervening mag- 
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netic fields. The magnetic field structure (both extragalactic 
and Galactic), along with the UHECR source distribution, the 
nature of sources (transient vs. steady), the energy spectrum 
at the injection, and the CR composition, are all "known un- 
known" factors that affect the distribution of the observed ar- 
rival directions. 

Though quite an extensive literature on simulation of 
UHECR propagation in extragalactic magnetic fields exists, 
very little is known about the strength and configuration of 
such fields. So far, direct evidence for the presence of extra- 
galactic magnetic fie lds has been found on ly in galaxy clus- 
ters (for a review see lCarilli & Taylorl2002l) . Faraday rotation 
measurements provide evidence for intracluster core fields in 
the range of 1 — 10 /iG. Outside clusters only upper limits at 
the level 1 — 10 nG are available. Extragalactic magnetic fields 
are ad hoc assumed to have a domain structure with a Kol- 
mogorov power spectrum and a uniform correlation length. 
Toy models assume a field strength ~1 nG in voids with a 
somewhat larger fi eld ~10 nG at the supergalactic plane (e.g., 
IStanev et alJl2003h . 

Recent simulations of magnetic fields in the intergalactic 
medium are more sophisticated. They take into account the 
growth of the magnetic fields from seed fields such that the 
resulting field strength traces the baryon density as the large- 
scale structure evolves. A more realistic extragalactic mag- 
netic field of this nature may result in significantly larger 
deflections than is expec t ed fro m a purely random field. A 
simulation by Sigl et al. (2004) used the Biermann battery 
mechanism to generate seed fields which were evolved, and 
then rescaled so that the magnetic field in the core of a simu- 
lated Coma-like galaxy cluster is comparable to the /iG fields 
as indicated by Faraday rotation measures. Simulations of 
large-scale structure formation and the build-up of magnetic 
fie lds in the intergala ctic medium have also been performed 
by Dolag et al. (2005). The basic assumption is that cosmo- 
logical magnetic fields grow in a magnetohydrodynamic am- 
plification process driven by the formation of structure from a 
magnetic seed field present at high redshift. The initial density 
fluctuations were constructed from the IRAS 1.2-Jy galaxy 
survey by first smoothing the observed galaxy density field 
on a scale of 7 Mpc, evolving it linearly back in time and 
then using it as a Gaussian constraint f or an otherwise ran - 
dom realization of a ACDM cosmology dMathis et alJl2002l) . 
As a result, the positions and masses of prominent galaxy 
clusters coin cide closely w i th thei r real counterparts in the lo- 
cal universe. iTakami et al.l d2006l) have used a magnetic field 
strength scaled with the matter density \B\ cx p 2 / 3 , where the 
distribution of galaxies is constructed using the IRAS PSCz 
catalog. The correlation length is taken to be 1 Mpc and the 
magnetic field is assumed to be represented as a Gaussian ran- 
dom field with a Kolmogorov power spectrum in each cube. 
The field is further renormalized to obtain ^0.4 //G in a cube 
that contains the center of the Virgo Cluster. 

The average deflection angle in a random field is (6d) ~ 
2.5° Z E^q -B_g(-Dioo^i) 1/2 , where Z is the particle charge, 
Dioo is the distance in units 100 Mpc, £?_g is the r.m.s. 
field strength in nG, E20 is the particle energy in units 
10 20 eV, and h is the correlat ion length in Mpc (e.g., 
IWaxman & Mira lda-Escude 1996). The average time delay 
correspond ing to the average defle ction angle is (r) ~ 
(9j)D/4c (lAlcock & Hatchettl I1978L their eq. [29]), where 
(O 2 ) — A{9 d) 2 /tt can be derived from 6 2 probability distri- 
bution (their eq. [23]), and c is the speed of light. This yields 



(r) ~ {6 d ) 2 D/nc - 2 x 10 5 Z 2 E w 2 B 2 _ 9 D'f 00 h yr. For a 
10 20 eV proton injected at D ~ 75 Mpc for characteristic val- 
ues of £_ 9 - 1 and h ~ 1, (6 d ) ~ 2° and (r) - 1.6 x 10 5 
yr which is comparable to the crossing time of the Galaxy and 
the characteristic timescale of jet lifetimes in Seyfert galaxies 
as discussed earlier, but is negligible compared to the galac- 
tic evolution timescale. Therefore, observations at different 
wavelengths show us nearly a snapshot of the sources at the 
time when the highest energy CRs were emitted. If a de- 
tected CR particle has been accelerated by a pc-scale jet, the 
jet will expand during the time delay to become larger, ^10 5 
lt-yr, and this may be seen as a more extended structure in 
the radio. The association of UHECR accelerators must cor- 
respondingly take into account such time delays and source 
evolution since observed photon signals come from later times 
than the epoch of UHECR escape from the source. 

The Galactic magnetic field is known much better than the 
extragalactic one. It can be determined from pulsar rotation 
and dispersion measures combi ned with a model for th e dis- 
tribution of free electrons (e.g.. lCordes & La zio 2003a|3). A 
large-scale field of a few /iG aligned with the spiral ar ms ex- 
ists, b ut there is no general agreement on the details dBeckl 
2001). Recent studies give a bisymmetric model for the 
large-scale Galactic magnetic field with reversals on arm- 
intera rm boundaries ( Bro wn et al.ll2.007h Han 200^ lHan et all 
2006). Independent estimates of the strength and distribu- 
tion of the field can be made by simultaneous analysis of 
radio synchrotron, CR, and 7-ray data, and these confirm 
a value of a few / iG, increasing towards the inner Galaxy 
(IStrong et al.l 12000. and Strong et al. in prep.). The mag- 
netic field in the halo is less known. Observations of the ro- 
tation measure of extragalactic radio sources reveal azimuth 
magnetic fields in the halo with reversed directions below and 
above the plane consistent w ith AO symmetry type (lHanl2008l : 
IWielebinski & Krausdll993l) . 

Because of their large Larmor radii >1 kpc (rx ~ 
10 5 i?2o/5_9 kpc), UHECRs propagating in the Galactic 
magnetic field are sensitive to the global topology of the 
field. The influence of the geometry of the Galactic magnetic 
fiel d has been studied in various source distribution scenar - 
ios (lAlvarez-Mufiizll2002l: Istanevlll997t ITakami et ai1l2006h . 
For a ~1 /iG magnetic field, the distance D ~ 100 kpc, and 
a correlation length I ~ 1 kpc, the average deflection angle 
is (Od) ~ 2.5°-E^ ) 1 , but the actual value depends on the ar- 
rival direction of a CR particle. Cen A is only ^50° away in 
longitude from the Galactic Center, and only ^20° from the 
Galactic plane, while Cen B is very close ~1° to the Galac- 
tic equator. Cosmic rays coming from either of these objects 
could be influenced by the stronger magnetic field near the 
Galactic plane (a few /iG vs. ~1 /iG in the Galactic halo) over 
tens of kpcs of their trajectory. This would provide a greater 
deflection than the relatively longer path length through the 
weak extragalactic magnetic field. Therefore, an association 
of Cen A and Cen B with more events in this region is plausi- 
ble. 

The UHECR source distribution is usually assumed homo- 
geneous or to follow the baryon density distribution. The for- 
mer case is relevant for energies below the photopion produc- 
tion threshold for proton injection where the energy losses are 
small and particles may come from cosmological distances. 
Since only a small fraction of the sky is covered with an extra- 
galactic magnetic field capable of deflecti ng UHECR particles 
by a significant angle (Dolag et al. 2005), the resulting distri- 
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bution of arrival directions is close to isotropic (Taka mTet alj 
2006&. If the observed energy of CRs is near the GZK cut- 
off, the sources are likely local. In this case the distribution 
of sources traces the baryon density distribution in the local 
Universe and the effective field acting on UHECRs should 
be considerably stronger, leading to larger deflection angles. 
The distribution of the deflec tion angles depe nds on the de- 
tails of the simulations, e.g.. ISigl et al.l (2004) pr edict s large 
deflection angles, -20° at 10 20 eV, while iDolag et al.l (2005) 
gives much smaller angles, but allows for angles >3° in a 
small fraction of the sky <0.01. Therefor e, the arrival di- 
rections of UHECRs (Abraham et al. 2007b) should correlate 
with the distribution of large deflections on a deflection map; 
such a correlation can be seen even from a by-eye comp arison 
with the deflection map given by Takami et al. (20061 their 
Figure 5). The sources of UHECRs should also be capable 
of producing lower energy CRs and 7-rays (and neutrinos) 
and, therefore, may be observed with the next generation 7- 
ray telescopes. 

Most of our discussion has described the situation if the 
UHECRs particles are protons. This is complicated further if 
the injected particles are CR nuclei since the deflection an- 
gles can be larger for a given magnetic field and the nuclei 
undergo photodisintegration processes on the CMB and ex- 
tragalactic infrared background frag menting into lighter nu- 
clei (e.g.. IStecker & Salamon1ll999l) . The UHECR chemi- 
cal composition is unknown and subject to consi derable de- 
bate. Results from the surface arrays AG ASA dShinozakil 
2006) and Yakutsk (Knurenko et al. 2008) and fluorescence 
detectors (e.g., Sokolsky & Thomson 2007, and references 
therein) indicate a trend toward proton dominated composi- 
tion at the highest energies. However, the Auge r Collabora- 
tion h as presented a fit to the elongation rate 14 (Ung er et alj 
120071) showing a heavier or mixed composition at the highest 
energies. These interpretations are complicated by the nec- 
essary reliance on hadronic interaction models which have 
to extrapolate cross section information beyond current ac- 
celerator energies, and indeed even details of the UHECR 
sources themselves can introduce degeneracy in the inter- 
pretati on of the data wi th different chemical compositions 
(e.g..lArisaka et alj|2007l). It has been argued r ecently (e.g., 
lDermeril2008tlFargionl2008tlHooper et alj|2008l) that the cur- 
rent anisotropy results can be explained if the composition has 
a significant component of light nuclei, 4 < A < 14, but this 
remains to be tested by further data. 

5. CONCLUSION 

A transition from an isotropic distri bution of ar rival di- 
rections of CRs above ~1 EeV dWatsonl 120081) to an 
anisot ropic distribution of th e highest-energy CRs above 57 
EeV dAbraham et alj l2~007bl) observed by Auger implies a 
change in the propagation mode of UHECRs in intergalactic 
and/or Galactic space. The ass ociation of the observed events 
with the supergalactic plane (Abraham et al. 2007b; Stanev 
2008) points to the sources tracing the supergalactic plane 
and matter distribution which correlates with AGN. However, 
as we have shown, almost all nearby (c?l < 150 Mpc) ac- 
tive galaxies found within the search radii of 3.2° around the 
UHECR events detected by Auger are typical for the local 
low-luminosity AGN of the Seyfert/LINER type. They are 
characterized by low-power and short jet activity, which is 

14 The elongation rate is the slope dXm^/d log E, where Jf max is the 
depth of shower maximum. 



substantially different from that observed in radio galaxies 
and quasars (as typically considered in the scenarios for accel- 
eration of UHECRs). Moreover, such selected low-luminosity 
AGN are expected to be quite common in the local Universe, 
with the estimated surface density 1.4 x 10 3 per steradian, 
when limited to a redshift z < 0.037. If the acceleration of 
UHECRs is indeed associated with jet activity, which is most 
likely, we conclude that the correlation with particular AGN 
is a coincidenc e. To distinguish betwe en the persistent and 
episodic (e.g., iFarrar & Gruzinovl l2008S) models of UHECR 
sources, future, more extensive analyzes have to take into ac- 
count details of the AGN radio morphology and spectral prop- 
erties, and may yield a correlation with a larger deflection an- 
gle and/or more distant sources. 

We emphasize that there is no complete all-sky catalog of 
nearby AGN. In addition, many "regular" galaxies when stud- 
ied at sufficient spatial resolution at different wavelengths 
show some (typically weak) level of the AGN-like activity. 
Hence, the confusion in classification of such sources in the 
literature, different databases, and catalogs. Thus, investi- 
gating the correlation of UHECRs with AGN based on some 
given particular AGN catalog may be tricky and even mean- 
ingless. In particular, usi ng catalogs of X-ray selected local 
active galaxies (as in, e.g., George et al. 2008) may give mis- 
leading results since the X-ray emission of Seyferts and LIN- 
ERs is produced by the accretion disks and the disk coronae, 
and therefore represents the accretion power of the active nu- 
cleus rather than the power of its jet. Although the most recent 
studies have found a correlation between the disk luminos- 
ity and the radio power of the unresolved nu cleus dHoll2002l ; 
iHo & PengfeOOU iPanessa et ai]l2006l 120071) . the former has 
no direct relation with the large-scale radio structures which 
are supposed to be capable of accelerating CRs up to the high- 
est energies. Besides, the present time X-ray luminosity of 
the disks may have nothing to do with the observed UHECR 
events because of the considerable time delay between the ar- 
rivals of particles and photons (see §4j; on the other hand, 
past UHECR acceleration activity that produced the observed 
UHECR events, if not episodic, has to manifest itself by ex- 
tended jets that we should be able to see now. As argued in 
this paper, the spectral and morphological properties of the 
jetted AGN which are selected as likely counterparts of the 
detected UHECR events should be considered in detail and 
compared with the properties of the parent population. 

Other possibilities include a few close sources with ex- 
tended jet/lobe structures, such as Cen A and Cen B, and rel- 
atively large deflections due to either stronger magnetic fields 
or due to the presence of heavy nuclei in the flux, or more 
distant sources. 

Observations with 7-ray telescopes, such as Fermi/LAT, 
HESS, MAGIC, and VERITAS may point to the class of 
sources able to accelerate particles to TeV energies, and are 
therefore potentially capable of accelerating particles up to 
EeV energies. Such sources could also produce TeV and 
UHE neutrinos. Taking into account the delay between the 
arrival times of 7-rays (neutrinos) and UHECRs, such obser- 
vations have to be interpreted with care: UHECRs may come 
from sources which are not generating TeV 7-rays anymore, 
or UHECRs that are accelerated in present day 7-ray emitters 
have not had time to propagate to us, yet. Meanwhile, Cen A 
and Cen B are two powerful nearby radio galaxies and, if they 
are indeed UHECR sources, 7-ray observations can provide a 
"current" picture at the time when the CRs were emitted since 
the overall time delay from propagation is very short. More- 
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over, Cen A is large eno ugh to be resolved by 7-ray instru- 
ments (e.g., Fermi/LAT, lAtwood et ai1l2008t iMcEnery efall 
2004). Therefore, observations with 7-ray telescopes may 
provide additional clues to the origin of UHECRs. 
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